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Abstract 

The Magnus series is an infinite series which arises in the study of 
linear ordinary differential equations. If the series converges, then the 
matrix exponential of the sum equals the fundamental solution of the 
differential equation. The question considered in this paper is: When does 
the series converge? The main result establishes a sufficient condition for 
convergence, which improves on several earlier results. 



1 Introduction 

The Magnus series is an infinite series which arises in the study of linear ordinary 
differential equations of the form y' — A{t) y, y{0) = yo, where y(t) is a vector 
and A{t) is a matrix. We assume throughout this paper that A{t) is a real-valued 
matrix, even though the Magnus series is valid in more general settings. 
The fundamental solution is defined by 

Y'^A{t)Y, Y{0)=1. (1) 

If A{t) is constant, then the solution of ([T]) is given by the matrix exponential 
Y{t) — e^*. This suggests the ansatz Y{t) — e^^^\ where f2 is a matrix function 
to be determined, for the nonautonomous equation ([T]). It turns out that f2(t) 
satisfies the differential equation 

oo 

n' = dexp^i {A{t)) = E if adf, {A{t)) , (2) 

fc=0 

where Bk denote the Bernoulli numbers (Bq = I, Bi = — i, B2 = ^, B^ = 0, 
etc.) and ada is the adjoint operator, defined by 

adniA) = [fl, A] = nA- Afl. 
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Magnus [H] applied Picard iteration on ^ to find an infinite series for f2: 



f A{T)dT-\[ f ' [A{T2),Ain)]dT2dn 
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A{T,),[A{T2).A{r,)] 



dr-i dT2 dri 



This series has since come to be called the Magnus series. We can write it as 

n{t) = ni{t) + n2{t) + n-sit) + n^it) + ■■■ , 

where the term r2„(i) is a sum of rt-fold integrals of rt — 1 nested commutators. 
Explicit expressions for the ^n{t) are given by Bialynicki-Birula, Miclnik and 
Plabahski [2] , Chacon and Fomenko [4] , and Iserles and N0rsett [9] . 

The Magnus series can be used to derive the Baker-Campbell-Hausdorff 
(BCH) formula for the product of two matrix exponentials. This formula states 
that e^^e^^ = e^ with 

B = A, + A2 + ^[Ai,A2] + ^[Ai,[AuA2]]+^[A2,[A2,A,]]+-- - . (4) 

Indeed, if we define the function 74( ■ ) by A{t) = A2 for t e (0, 1) and A{t) = Ai 
for t G (1,2), then the solution of (jT]) at time i = 2 is e^le^^ and the BCH 
formula ([4]) is the Magnus series ([3]) for this particular choice of ^( • ). 

Magnus derived the expansion fS]) in the context of quantum mechanics, 
where A{t) is a skew-Hermitian matrix. Because the Magnus series is contructed 
from commutators, i}{t) is also skew-Hermitian and e^*-*' is unitary, just like the 
fundamental solution of the original differential equation ([T]). More generally, if 
A(t) is in some Lie algebra, then D,{t) will be in the same Lie algebra and e^*^*^ 
will be in the corresponding Lie group. 

In the 1990s, Arieh Iserles and Syvert N0rsett were among the group of 
mathematicians that established the discipline of Geometric Integration. This 
field concerns itself with numerical integrators that respect the geometric struc- 
ture of differential equations (see, e.g., Hairer, Lubich and Wanner [7j for an 
introduction). Iserles and N0rsett were looking for a way to integrate ^ such 
that the numerical solution evolves on the Lie group if the matrix A{t) is in the 
Lie algebra. Unaware of Magnus' work, they rederived the Magnus series. They 
realized that after truncating the infinite series and approximating the integrals 
by quadrature, a method arises that respects the Lie-algebraic structure of the 
equation (see Iserles and N0rsett [9] and Iserles for details). 

Since the Magnus expansion is an infinite series, it is natural to ask whether 
it converges. Indeed, Magnus himself gave an example in which the series does 
not converge [12]. The question of convergence is the subject of this paper. The 
main result (Theorem [3]) gives a sufficient condition for convergence. 
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2 Previous results on the convergence 



Due to the complexity of the expansion Q, several proof strategies have been 
used to derive bounds on the terms which in turn have led to many different 
convergence estimates. Magnus |12) gave no convergence estimate but stated 
that for sufficiently small t the series converges. By the known bound 



one easily arrives at the conclusion that the Magnus series converges whenever 



However there are several improvements on this bound. In the following, r 
denotes a number for which the following statement holds: 



In the field of quantum physics there has been some interest in the con- 
vergence issue. In 1966, Pechukas and Light [T7] consider particular quantum 
systems and find convergence conditions, althought these are not of the general 
form we are considering here (see also Fernandez [5], Klarsfeld and Oteo [TT], 
and Salzman [18j). 

In 1976, Karasev and Mosolova [10] cite a bound r = i log 20 Agrachev 
and Gamkrelidze, working in the field of control theory, mention in their 1981 
paper [1] a result by Vakhrameev stating that r = 1.086880 In 1998, this 
bound was rediscovered independently by Blanes, Casas, Oteo and Ros [3] and 
by Moan [T3], using different methods. Vela [52] states in 2003 that this is a 
sharp result. 

In 1991, Chacon and Fomenko [3] found an alternative expression for f7„. 
They used this expression to prove that r = 0.57745. 

A few years earlier, in 1987, Strichartz [12] had rediscovered the explicit 
expression for f2„ found by Bialynicki-Birula et al. [5], stated in terms of Lie 
brackets. He used this to prove r = 1. The same result was found independently 
by Vinokurov [23] in 1997. Finally, Moan and Oteo [15] derived the bound r = 2 
(the best result at the moment) by similar techniques, except that they avoided 
the use of commutators as they seemed to introduce unneccesary complications 
in the convergence bound. 

Moan [13] found a condition for existence of a real logarithm for real A with 
r = TT. It was however unclear if this condition is sufficient for convergence of 
the series expansion. Theorem [3] answers this question affirmatively. 

^V.S. Varadarajan |211 p. 119] gives implicitly the same result for the BCH formula. 
^In 1977, Suzuki gO] derives this bound for BCH. 

•^Newman et al. I16| establish the same bound for the BCH formula in 1988. 




/o*P(r)||2dr<log2Q 



If / ||^(r)||2 dr < r then the Magnus expansion converges. 

"'0 
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3 The existence of a real logarithm 



The fundamental solution Y{t) of the differential equation ([T]) is an invertible 
matrix. A theorem by Gantmacher [6] states that every invertible matrix has a 
logarithm, and hence there exists an n{t) such that e^'*) — Y{t). However, the 
logarithm may fail to be real, even if Y(t) is real. For example, the matrix 

6eM\{0}, (5) 

does not have a real logarithm. Since all the terms in the Magnus expansion ([3]) 
are real if the original differential equation ([1]) is real, we conclude that the infi- 
nite series cannot converge to a logarithm of Y{t) if Y{t) has no real logarithm. 
The question we are studying in this section is therefore: Does Y{t) have a real 
logarithm? 

The following lemma is easily proved, for instance by factoring $ — VJV^^ 
where J is in Jordan form. 

Lemma 1. Suppose that the invertible matrix $ has no negative eigenvalues, 
that is, suppose that the eigenvalues of ^ are contained m C \ (— cx),0]. Then 

poo -, 

log$ = (ci>-I) / {^l + q>)~Ufi. (6) 

Jo -L + M 

It follows that if Y{t) is real and has no negative eigenvalues, then the logarithm 
of Y{t) is real as well. The next result (which can be found in [T3]) gives an 
easy condition under which Y{t) has no negative eigenvalues and hence a real 
logarithm. 

Theorem 2. Let A{t) be a real integrable matrix, and let Y{t) denote the so- 
lution ofY' = A{t)Y, y(0) = I. // J* \\A{T))\\2dT < n, then Y{t) has a real 
logarithm. 

Proof. Choose an arbitrary vector yo, and consider the vector y{t) — Y{t)yo 
satisfying y' = A{t)y and y(0) — j/o- Let y{t) = y{t)/\\y{t)\\2 denote the unit 
vector in the direction y{t). Then 

y'{t)^{iut)h)m + \\ym2y'it). 

Since y{t) and y'{t) are orthogonal, it follows that 
Therefore, we have 

\Wm2\\ym2<\\y'it)h<\\Ait)h\\yit)U 

and thus ||y'(f)||2 < ||^(i)||2- Integrating this inequality, we find the bound 

'\\y'ir)hdr< f \\A{r)hdr. 
Jo 



-1 

b -1 
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The left-hand side is the length of the curve swept out by the unit vector ^'(t) 
when T G [0,i\. Therefore, the angle between y{t) and y(0) is smaller than tt 
if /p \\A{T))\\2dT < TT. Since y{t) = Y{t)y{0), this implies that Y(t) has no 
negative eigenvalues, and hence it has a real logarithm. □ 

Example[2]in the forthcoming Section [5] shows that the constant tt in Theorem[2] 
is sharp. 

4 Proof of convergence 

Theorem [5] gives a condition for Y{t) to have a real logarithm. The following 
theorem states that under the same condition, the Magnus series converges to 
this logarithm. 

Theorem 3. Let A{t) be a real integrable matrix, and let Y{t) denote the solu- 
tion ofY' = A{t) Y, Y{0) = I. // /q ||A(t)||2 dr < tt, then the Magnus series © 
converges and its sum Q,{t) satisfies e^^*-*' = Y{t). 

Proof. We write the Magnus series as VL{t) = Y^'^=i^n{i)^ where the term 
f2„(t) is a sum of n-fold integrals of n — 1 nested commutators. If we now 
introduce a new parameter k, and we replace A{t) by K.A{t), then the Magnus 
series becomes Q,{t] k) = The idea is to fix t and consider the 

function / defined by /(k) = logY{t;K) where Y{t;K) denotes the solution of 
Y' = fiA{t)Y, Y{0) = I with n e C. We will show that / is analytic and that 
the Magnus series is the Taylor series of this function around k = 0. We can 
then use the standard result from the theory of complex functions which states 
that the Taylor series of a function converges in a disc if the function is analytic 
in that disc. 

Set 7 = /o ||^(''')||2 rf'''- As stated in Section [2l it is easy to show that if 
|«;| < i log 2, the Magnus series converges and its sum Q{t; n) satisfies e^^'^*'''^ = 
Y{t;K). Hence, the power series V,{t; k) coincides with /(k) for |k| < ilog2, 
and the Magnus series is the Taylor series expansion of / around k = 0. 

We say that a matrix-valued function (like /) is analytic if all the matrix 
entries are analytic functions. We now want to prove that / is analytic in the 
disc with radius ^. Firstly, the fundamental matrix Y(t]K) is analytic as a 
function of k. The proof of Theorem [5] shows that Y{t; k) has no eigenvalues in 
(—00,0] if \k\ < ^, so the logarithm is given by ©. Hence, the derivative of / 
is given by 

= j^Yit; k) / — - (mI + y (t; k)) dfi 

OK Jq f-L -\- 1 

~Yit;K) / {^l + Y{t■,K)y^—Y{t■,K){^ll + Yit;K)y'd^l. 

Jo + -L OK 

The right-hand side is well-defined, thus proving that the function / is analytic 
in the disc with radius -. Hence, we can expand the entries of the matrix /(k) 
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in a power series, and this series will converge provided that |k| 
power series is precisely the Magnus series. 



< ^. But this 
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5 Examples 

In this section, we study some examples to investigate the connections between 
the condition \\A{T)\\2dT < tt, the eigenvalues of the fundamental solution, 
the existence of a real logarithm, and the convergence of the Magnus series. 

Example 1. The following simple example suffices to show that the condition 
||A(t)||2 dr < vr is not a necessary condition. Consider the equation Y' = AY 
where A = [ _?i J ] is a constant matrix. The fundamental matrix is 



Y{t) 



cost 
sint 



— sint 
cos t 



When t = TT, this matrix has a double eigenvalue at —1. Furthermore, we have 
Jq \\A{t)\\2dt ~ TT. Nevertheless, Y{t) has a real logarithm for all t, including 
t — TT, because Y{t) = e*^. The Magnus series converges for all t: the terms in 
the series are fti — A and Ofc = for A: > 1. The critical point here seems to 
be that the double eigenvalue at —1 is not defective (meaning that its algebraic 
multiplicity equals its geometric multiplicity). 

Example 2. This example, taken from Moan [14], shows that the condition 
/p ||A(r)||2 dr < TT is sharp, in the sense that the constant it on the right-hand 
side cannot be replaced by a bigger constant. Consider the equation Y' = A{t)Y 
with 



Ait) ^ - 



sin 2t 
1 — cos 2t 



A simple computation shows that \\A{t) 
the differential equation is given by 



Y{t) 



t sin t + cos t 
sin t — t cos t 



L — cos 2t 
— sin 2t 

dr = vr, and that the solution of 



— sint 
cost 



Hence l"(7r) — [ , which is of the form ([5]), and therefore F(7r) does not 

have a real logarithm. This implies that the Magnus series diverges for t — tt. 

Example 3. This example shows that the Magnus series may diverge even 
though the fimdamental matrix Y{t) has a real logarithm. Take 



Ait)^ 



t 

-1 



The Magnus series is 

r2(t) 



r2t it^i 


+ 





-it^l 

4'- 


+ 





0' 


+ 


^ so'' 


-t 
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On the other hand, the solution of Y' = A{t)Y is 



Y{t) 



1 2t 
9^^ 







The logarithm of this solution is 



\ogY{t) 



2t fit) 
-t 



where /(t) — 



ie^* - (t + 3i2)e" 



3(e- 



e*) 



So, the Magnus series is the Taylor series expansion of log Y{t) around t = 
(this is not true in general). The function / has a pole at t = |7ri, thus the 
Magnus series converges up to < = |7r. 

However, Y{^n) has eigenvalues at e'^^/^ and e^'^^'^, so Yi^n) has a real 
logarithm. Nevertheless, the Magnus series diverges. The reason is probably 
as follows. As in the proof of Theorem [31 let Y{t;K) denote the fundamental 
solution of Y' — KA{t)Y. The eigenvalues of Y{t; k) are e*^' and e"^''*. When 
K — i, these eigenvalues move in a circle around the origin and collide when 
t = ^TT. This collision causes the Magnus series to diverge at i = Itt. 



Example 4. The preceding examples all have some structure which allowed us 
to determine where the Magnus series starts to diverge. We close this section 
with a more-or-less randomly chosen example. 

Consider the equation Y' = A{t)Y, where the matrix A is defined by 



A{t) 



-t 3t 

-t -3 + 2t + 3 

3 -2t 

~t + 3 2t^ -it -it - 2 



-M^ + t + 2 


-<2-3 
-t + 2 



(7) 



Since A is polynomial in i, all the commutators in the Magnus series are also 
polynomials and hence the integrals can be computed exactly. The first two 
terms of the series are given by 



Qi = 



^2 = 



1+3 



_ 1+2 



1+3 



it 

\t^ + it 



60 ^2 ^4 
60 ' 4 2 



3+2 

2"- 
-it 



2+3 



3 .2 





\t^ +t^ +it 

- t^ 



42 





-t^ + + 2t 



3^2 



-t='^^t +-t 



1*3 



2t 



it 

t^ + 2t 



10 ^ 2 ' 12 ^ 

J-fS , i(4_ 3 ,3 
1 (5 , 2j4 5.(3 



_ J_tS+i fi 

60 ^3 
.J_t5 , 3j4 , 1(3 
12 ^4 ^3 



We computed the first thirty terms of the Magnus series with the help of a 
computer algebra system. The recursive formulas given by Blanes, Casas, Oteo 
and Ros [3j proved to be useful for this purpose; other formulations require a 
very long time to evaluate. In Figure [Tl we plot the sums of the first fifteen, 
twenty, twenty-five and thirty terms of the Magnus series ([3]) . It seems that the 
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Figure 1: The left plot shows the partial sums of the first 15, 20, 25 and 30 
terms of the (1,1) element in the Magnus expansion ([3]), when A is as given 
in ([T]). The dash line indicates t = 0.733. The right panel shows the same for 
the (2,3) element. 

series starts to diverge between t — 0.7 and t — 0.8, though it is of course not 
possible to pinpoint the location precisely. 

Figure [2] shows how the eigenvalues of the fundamental solution Y(t) move 
around as t increases from to 1. The eigenvalues do not collide in the left 
half-plane, so Y{t) has a real logarithm for t G [0,1]. However, as the right 
panel in the figure shows, the eigenvalues of Y{t; e*"*) with a* = 1.805 ... do 
coUide. The collision takes place at = 0.733 . . . and A* w -0.485 + 0.0249i. 
The value of is shown by the dash line in Figure [TJ It seems plausible that 
the Magnus series starts to diverge around this point. 

Incidentally, the condition ||A(r)||2rfT < tt is satisfied for t < 0.56.... 
However, the Magnus series continues to converge for slightly larger values of t, 
showing again that this condition is not necessary for convergence. At t = t^,, 
the integral is approximately 4.36. 

6 Conclusion 

The main result of the paper is Theorem [31 which states that the Magnus series 
converges if ||^(r)||2 c?t < tt. This condition is in the same form as various 
earlier results mentioned in Section [2l Example [2] shows that our result is sharp 
in the sense that the constant tt is the largest number for which the result 
holds. However, the other examples show that the condition is not necessary 
for convergence. 

The proof of Theorem [3] shows that the Magnus series can be considered as 
the Taylor series oilogY(t; k) around k = 0, and hence the radius of convergence 
is determined by the nearest singularity. Lemma [1] implies that there are no 
singularities if the eigenvalues of Y(t; k), which start at 1 when k = 0, do not 
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1 



0.2 




Figure 2: The left plot shows the location in the complex plane of the eigenvalues 
of the fundamental solution of Y' = A{t)Y with A as in ([ll), for t G [0, 1]. All 
four eigenvalues are at 1 when t = 0. As t increases, two eigenvalues move to 
the right, while the other two form a complex pair. The small circles show the 
locations at t = 0, 0.1, 0.2, . . . , 1. It is clear that no eigenvalues collide on the 
negative real axis for t G [0, 1]. 

The right plot shows the eigenvalues for ¥' = e'^"' A{t)Y with a, « 1.805. Two 
eigenvalues encircle the origin and collide near the arrow. 

cross the negative real axis as k moves in the unit disc. However, the choice of 
the negative real axis is arbitary; a formula similar to ^ holds for any other 
branch cut of the logarithm. 

Examples [3] and m suggest that divergence of the Magnus series is associated 
with eigenvalues of Y{t] n) encircling the origin and colliding as k moves from 0. 
In this situation, the eigenvalues are on different sheets of the Riemann surface 
of the logarithm when they collide. The authors feel that this conflict leads to 
the divergence of the Magnus series. However, as example [1] shows, not every 
collision of eigenvalues leads to divergence of the Magnus series. If the multiple 
eigenvalue at the collision is not defective, then the eigenvalues retain their 
identity throughout the process and no conflict arises. 

Of course, this is not a proof, but we think that it may be possible to make 
it rigorous using the correct formalism. We are thus led to propose the following 
conjecture, which gives a necessary and sufficient condition for convergence. 

Conjecture 4. Let Y{t; k) denote the solution of Y' ~ KA{t)Y , Y{0) — I. 
Denote the eigenvalues ofY{t;K) by Xn{t;K), where the eigenvalues are to be 
ordered so that A„ is a continuous function oft. Let i» be the smallest t > Q for 
which there exists a k G C with \k\ — 1 such that there is a multiple eigenvalue, 
say Xi(t, k) — Xj{t, k) with i ^ j , for which the geometric multiplicity is smaller 
than the algebraic multiplicity, and the loop 

{A,;(t,k) I r G [0,i]}U{Aj(T,K) I T G [0,i]} 

encircles the origin. Then the Magnus series converges if and only if t < t^. 
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This conjecture suggests two tasks. The first is obviously whether we can find 
a proof for this conjecture. However, even if the conjecture is true, it is difficult 
to apply in practice because the condition is not easy to check. Therefore, the 
second task is to find a more practical condition for convergence. 
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by EPSRC First Grant GR/S22134/01. 
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